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ABSTRACT 

Using semi-analytical models of galaxy formation, the origin of boxy and disky ellip- 
tical galaxies is investigated. We find that the simple scenario, motivated by N-body 
simulations, in which the isophotal shape is only dependent on the mass ratio of the 
last major merger, is not able to reproduce the observation that the fraction of boxy 
and disky ellipticals depends on galaxy luminosity. The observations can however be 
reproduced with the following reasonable assumptions: (i) equal-mass mergers lead 
to boxy ellipticals and unequal-mass mergers produce disky ellipticals (as motivated 
by N-body simulations) (ii) major mergers between bulge-dominated galaxies result 
always in boxy ellipticals, independent of the mass ratio, (iii) merger remnants that 
subsequently accrete gas leading to a secondary stellar disk with more than 20% the 
total stellar fraction are always disky. This scenario indicates that the isophotal shapes 
of merger remnants are sensitive to the morphology of their progenitors and subsequent 
gas infall. Boxy and disky ellipticals can be divided into two subclasses, depending on 
their formation history. Boxy ellipticals are either formed by equal mass mergers of 
disk galaxies or by major mergers of early-type galaxies. We find that disky ellipticals 
are indicators of unequal-mass mergers or late gas infall. Disky ellipticals with high 
luminosities are preferentially 1-component systems that result from unequal mass 
mergers whereas low-luminosity disky ellipticals are more likely to harbour secondary 
disk components. In addition, the fraction of disky ellipticals with secular disk compo- 
nents should increase in regions with higher galaxy densities. Taking into account the 
conversion of cuspy cores into flat low-density cores by black hole merging we find that 
disky ellipticals should contain central density cusps whereas boxy ellipticals should 
in general be characterised by flat cores. Only rare low- luminosity boxy ellipticals, 
resulting from equal-mass mergers of disk galaxies could have power-law cores. 
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1 INTRODUCTION 

Elliptical galaxies are old stellar systems that are believed 
to have been formed by major mergers of disk galax- 
ies, preferentially at high red s hift ijXoomre fc Toomrelll972l : 
ISearle. Sargent, fc Baenuold 119731: ICimatti et alJ |2004) . 
This "merger hy pothesis" has been tested using numerical 
simulations (e.g. iGerhardllTQS l: Ncgropontc & White 1983; 
iBarnedligS 3; lHerrMuis^n992ri which demonstrated consis- 
tently that the global properties of equal mass merger rem- 
nants resemble those of ordinary slowly rotating massive el- 
liptical galaxies. 

More recently it has become clear that ellipticals are 
more complex than originally thought. Isophotal fine struc- 
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tures have bee n detected tha t correlate well with kinemati- 
cal properties jBendeilllQsl) . Faint ellipticals are isotropic, 
fast rotators with small minor axis rotation. They are called 
disky as a Fourier analysis of their isophotal deviations 
from perfect ellipses leads to positive values of the fourth 
order coefficient 04. Disky ellipticals might contain sec- 
ondary, faint disk components which contribute up to 30% 
to the total light in the galaxy, indicating dis k-to-bulge ra- 
tios t hat overlap with those of SO-galaxies jRix fc Whitd 
I1990I: IScorza fc Bendeil ll995!). In contrast to disky ellipti- 
cals, boxy systems are characterised by negative values of 04. 
They rotate slowly an d are supported by velocity anisotropy . 
Their cores are flat jLauer et alJ Il995l : iFaber et aP 119971) 
and reveal complex internal kinematics. The distinct physi- 
cal properties of disky and boxy elliptical galaxies indicate 
that both types of ellipticals experienced different forma- 
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tion histories. Interestingly, most of the massive ellipticals 
are boxy, while 2 /3rd of the lower -mass ellipticals are disky 
ijBender. Burstein. fc FaheAhm'j ) . 

In order to understand the origin of boxy and disky el- 
lipticals within the framework of the major merger scenario, 
numerical merger remnants of dissipationless merg ers have 
been a n alysed in detail e.g. bvlH ernquistj (|1993|'l.,Hevl et alJ 
iim4). iLiina-Neto fc CombesI (ll^95^ ■ iNaab et alJ t99^) 
and lBend^7^Bara^Til200ol) . In addition the role of gaseous 
dissipation on the the isophotal shape of remnants has been 
invest ig ated in various s tudies e.g. byjBarncs fc H ernauisd 
JlQQeTl . lBekki fc Shioval lll997D and ISoringel 1.200a) . These 
studies find in agreement with each other that gaseous 
dissipation leads to enhanced diskyness in the remnant. 
[Sckki fc Shiova ( 1997) simulated a number of mergers be- 
tween bulgeless disk galaxies with purely gaseous disks 
which depending on the rapidity of star formation tend 
to be more likely disky (boxy) for weak (strong) star for- 
mation efficiency. This result can be understood in terms 
of gas that has lost large parts of its angular momen- 
tum and settles down in the centre of the remnant, in- 
creasing the potential well and thereby destabilising box - 
orbits passing the central region jBarnes fc Hernauistll99'3) . 
The amount of gas able to settle down into the centre 
is strongly correlated with the efficiency by which gas is 
transformed into stars whic h explain s the trend observed 
in the simulations. However, ISpringeil ll2000ll demonstrated 
that this effect is mu ch weaker in the presence of a stel- 
lar bulge componen t . [Kauffmann fc Haehneld ||200(J) and 
iKhochfar fc Burkerd ^2003^ showed that the role of dissi- 
pation in the formation of massive spheroids is weaker com- 
pared to its role in the formation of low mass spheroids. Fur- 
thermore massive spheroids tend to form mainly by mergers 
of pro genitors already consisting of massive bulge compo- 
nents iKhochfar fc Burkerdl2003) . In the following, we will 
compare our results to the giant and parts of the inter- 
mediate elliptical sample with absolute B-band magnitude 
Mb ^ -18.5 of Bender, Burstein, & Fabcr (1992) allowing 
us to neglect the influence of dissipation during the merger 
process on the isophotal shape. 

Recentlv. iNaab fc Burkerd J2003D (for a summary see 
also lBurkert fc Naabll2003D completed a large survey of dis- 
sipationless merger simulations of disk galaxies, adopting 
a statistically unbiased sample of orbital initial conditions 
with mass ratios of 1:1 to 4:1. They showed that unequal 
mass 3:1 to 4:1 mergers lead to fast rotating disky ellipticals, 
in excellent agreement with observations. In contrast, slowly 
rotating, pressure supported ellipticals formed in equal mass 
1:1 to 2:1 mergers of disk galaxies while 2:1 to 3:1 mergers 
lead to remnants with mixed isophotal shapes. 

In this paper we investigate whether the scenario of 
disky and boxy elliptical formation based on mass ratios 
in major mergers is in agreement with the observed distri- 
bution of ellipticals with different isophotal shape. Using 
semi-analytical simulations we predict the ratio of boxy and 
disky ellipticals as a function of their luminosity. In section 
2 fc 3 we introduce the semi-analytic approach we use and 
the way in which we assign ispohotes to modelled ellipticals. 
Section 4 compares our model predictions with observations 
and in section 5 we introduce two new sub-classes of ellipti- 
cal galaxies followed by a concluding section. 



2 FORMATION OF ELLIPTICALS 

We follow the formation and evolution of elliptical galaxies 
in the context of semi-analytical modelling of galaxy forma- 
tion. The history of dark matter halos is traced using the 
merger tree proposed by Somorvillo fc Kolatt ( 1999) and the 
baryo nic physics is treated as described e.g. in lSoringel et alJ 
(|2001j) . The dark matter merger trees are resolved down to a 
minimum mass of Mmm ~ 10^" Mq and calculated for a set 
of different present day dark halo masses {Mo = 10", 10", 
lO^'' Mq). Larger halo masses correspond to more dense en- 
vironments, as e.g. Mo = lO" and Mo = lO" M© cor- 
responds to a present day cluster and field environment, re- 
spectively. The resolution limit we adopt does not infiuence 
our results as we will show in the next section. We assume 
a ACDM model with cosmological parameters Q,m ~ 0.3, 
Qa = 0.7, h = 0.7, Qm/^b = 0.16 and as = 0.9. 

The galaxy formation paradigm inherent to semi- 
analytic models assumes that the first galaxies to form are 
bulge less disk galaxies (see however iD'Onghia fc Burkerd 
l2004f) which form s tars in centrifugally s upported disks 
of cooled gas (e.g. IWhite fc FrenkI [19911) . Subsequently 
spheroids, elliptical galaxies and bulges, are formed in ma- 
jor mergers with mass ratios Mga i i/Mggi.i < 3.5 (Maai.i ^ 
Mga i 2) of galaxies as proposed bv[T oomre fc Toomr3 (^T^) 
(e.g. iKauffmann et alJll999l: ISpringel et alJl200lD. We adopt 
a sche me identical to the one used in e.g. ISoringel et al] 
i200lf) in which each galaxy in the semi-analytical simula- 
tion consists of two stellar components, a disk and a bulge 
component. 

2.1 Disk Component 

The disk component grows through hot gas that radiatively 
cools from the halo region and settles into the equatorial 
plane of a galaxy where it is transformed into a secular stel- 
lar disk. The co oling rate is c alculated using the prescrip- 
tion described in lSpringel et a l. (2001). Additionally we al- 
low disks to grow by accretion of the cold gas in the disks of 
satellite galaxies in minor mergers. We neglect any further 
possible ways of growing disks and comment on that in the 
next section. 

2.2 Bulge Component 

The bulge component grows through major mergers of 
galaxies. We calculate the timescale for galaxies to merge us- 
ing the approach in Kauffmami et al. ( 1999) with the mod- 
ification of using th e approximation of th e Coulomb loga- 
rithm as proposed in lSoringel et alJ (I2OOII) . Mergers disrupt 
the proge nitor disks as seen in various numerical simula- 
tions (e.g. lBarnes fc Hernauisdll992HBurkert fc Naabll2003l. 
and reference therein) and relax to a spheroidal distribution. 
During the merger any cold gas in the disk of the progeni- 
tor galaxies is assumed to be funnelled into the centre of the 
remnant where it ignites a star burst which transforms all of 
the cold gas into stars contributing to the spheroidal compo- 
nent (e.g. jCauffmarm ct al. J_9.99,; Spr ingel et al. 2001,, and 
reference therein) . The second assumption is certainly a sim- 
plification of what might happen since we neglect that not 
all of the cold gas is funnelled to the centre but some fraction 
of it can e.g. settle down in an extended disk which continues 



KrowinR inside out by fresh supply of gas from t idal tails ( e.g. 
Barnes fc Hernauisjlll99lt iMihos fc Hernauistl [l996: Barne3 
2OOII. l2002h . iBarnesI lj2002) results indicate that 40% - 80% 
of the initial gas could end up in the central region of the 
remnant and be consumed in a star burst. The exact amount 
is somewhat dependent on the merger geometry and on the 
mass ratio of the merger. Unfortunately, a large suryey in- 
vestigating the gas inflow to the centre of merger remnants 
is not available to date so that we use the simplified ap- 
proach. The prescription for the faith of the cold gas we 
adopt results in an overestimate of the spheroid masses and 
an underestimate of the secondary disk components in our 
model, which is not very significant for massive ellipticals 
since the y are mainly formed in relatively weak dissipative 
mergers Kauffmann fc Haehnehl2000l : iKhochfar fc Burkerii 
l2003r . Another simplifying assumption is that we neglect 
the feeding of super massive black holes in the centre of 
the remnant or feedb ack effects on the gas fr o m th e cen- 
tral source. However, iHaehnelt fc KauffmannI i2000h esti- 
mate that a cold gas mass fraction of less than 1% accreted 
onto the black hole is sufficient to recover the M. — a re- 
lation and we therefore neglect this effect. Furthermore we 
assume, that the bulge components of galaxies grow in minor 
mergers by the stars of infalling satellites. 
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Figure 1. Upper panel: Ratio of boxy to disky ellipticals for 
dilTerent magnitudes. Symbols refer to runs with different mass 
resolution. The labels in the graph indicate the halo mass at a 
redshift of z = (in this case Mo = M0) followed by the 

minimum mass (Mmin = 10^°, 5 X 10^, 10^, 5 X 10* Mq). Lower 
panel: Ratio of boxy to disky ellipticals for different magnitudes 
found in the simulation. Symbols refer to runs with different dark 
halo masses at 2; = 0. The labels in the graph indicate the halo 
mass at a redshift of 2 = (in this case Mq = lO^^, lO", lO^^ 
Mq) followed by the mass resolution {M,ni„ = 10^" Mq). 



2.3 Morphological Classification 

The above description for the growth of different stellar 
components allows for morphological transitions of individ- 
ual galaxies. All galaxies start off as bulgeless disk galaxies 
which later transform into elliptical galaxies in major merg- 
ers, and then possibly grow new secondary disks through 
the combined effect of cooling of hot gas from halo regions 
and cold gas accretion in minor mergers. In this way, pro- 
vided that enough cold gas is accreted, the ellipticals could 
also become bulges of early-type spirals until the next ma- 
jor merger which destroys the disk again and creates a new 
elliptical galaxy when the process of disk growing starts 
again. The knowledge of the bulge and disk components al- 
lows a morphological classification of modelled present day 
galaxies using the correlation between the Hubble type T 
and the _B— band bulge-to - disk r atio of galaxies presented in 
ISimien fc de Vaucouleur^ il986l) . Elliptical galaxies are clas- 
sified by T < -2.5 and have AM = Mb, bulge - MB,tot ^ 
0.55. Our definition of morphology bares some danger with 
respect to lenticular galaxies. Generally semi-analytic mod- 
els assume —2.5 <T < 0.92 for SO galaxies, identifying them 
as a transition type between ellipticals and spirals, in the 
process of regrowing a s econdary disk after a major merger 
ilKauffmann et alJl99gl) . Alternative formation mechanisms 
e.g. involve either the unequal mass mergers of two spiral 
galaxies (Bckki 1998; Barnes 1996; Bcndo & Barnes 200(J) 
in which the larger progenitor's disk basically survives the 
merger, or the growing of a secondary disk by material from 
tidal tails llMihojl20ol) . As a consequence, our population 
of ellipticals formed in mergers should be polluted with SO 
galaxies. At this stage, a clear separation between modelled 
disky ellipticals and SO galaxies is not possible and is hin- 
dered by observational ambiguities due to projection effects. 
We therefore do not distinguish between disky ellipticals and 
SO galaxies for T < -2.5. 



3 ASSIGNING ISOPHOTES 

Following the work of "Barnes' (l99^, iNaab et al.1 
and Naab & Burkort ( 2003) wc assign isophotes to ellip- 
tical galaxies according to the mass ratio of their last 
major merger. Last major mergers with mass ratio 1 ^ 

Mgal.l/Mgata < 2 and 2 < Mgal,l/Mgal.2 < 3.5 (Mgat,! 

^'Igai.2) lead to boxy and disky ellipticals, respectively. 

Fig.0shows the ratio of boxy to disky elliptical galaxies 
depending on their present day _B— band magnitude, assum- 
ing a final dark halo mass of 10^^ M©. The ratio is almost 
constant with Nboxy/Ndiaky ^ 0.5 over the presented mag- 
nitude range. This result is a consequence of the self-similar 
build up of elliptical galaxies by major mergers on all mass 
scales. We mentioned above that we follow dark matter his- 
tories back until progenitor halos drop below Mmin = 10^" 
Mq . The different symbols in the upper panel Fig. show 
the results of a resolution study of this minimum mass. By 
increasing the resolution (i.e. decreasing Mmin) we also in- 
crease the number of resolved small mass satellite galaxies 
in the overall galaxy population. As shown, the results do 
not change which is due to small mass satellite galaxies not 
significantly taking part in the last major major mergers at 
these luminosity scales. In the following we will use a mass 
resolution of Mmin ~ 10^" Mq which is a good compromise 
between efficiency and accuracy. 

It is well known that, observationally, the fraction of el- 
liptical galaxi es increases wi th the galaxy density of the envi- 
ronment (e.g. lDressleill980|). a trend also seen in simulations 
of galaxy formation (e.g. Ispringel et al ] |200ll and reference 
therein), which is attributed to the increased merger rate. 
The fraction of major mergers and its redshift evolution are 
very sensitive to environmental effects jKhochfar fc Burkerti 
l2001ri . However, one is not expecting to find the relative frac- 
tion of major mergers for a given mass ratio to change with 
the environment. The lower panel of Fig. 1 shows the de- 
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Figure 2. Ratio of boxy to disky ellipticals for different mag- 
nitudes found in the sim ulation. Filled circles show the data 
of iBendcr. B urstein. fc F aber 1 1992) and filled squares the re- 
sults for the run with Mq = 10^^ Mq and mass resolution 
Mmin = 10^° Mp). 



pendence of the ratio of boxy-to-disky ellipticals on the final 
dark halo mass which corresponds to different environments. 
Large final masses represent higher density environments, 
like clusters, low-mass represent a field environment. There 
is no significant environmental dependency, suggesting that 
the mix of boxy to disky ellipticals is universal. 



4 COMPARISON WITH OBSERVATIONS 

The next step is to compare the model predict ions to ob- 
servations by iBender. Burstein. fc Faberl (|l992h who pre- 
sented a first systematic study of the correlation of the 
isop hotal shape with o ther physical properti e s of e llipti- 
cals. iNaab et alJ and iNaab fc Burkerd J2003D used 

the sa me reduction routines as lBend^r^Burstein^g Faberl 
il992h in their merging simulations to derive the isophotal 
properties of their remnants. It is therefore safe to assume 
that we do not introduce any ambiguities in our definition 
of boxy and disky ellipticals by adopting their predicted de- 
pendence of isophotal shape on the mass ratio of the merger. 



4.1 A First Test 



The data of lBender. Burstein. fc Faberl lll992h shows a trend 
that the most massive ellipticals are mainly boxy, while the 
least massive ones are mainly disky. We compare this data 
with the model predictions in Fig.|5| This comparison shows 
a clear failure of the theoretical model as it does not repro- 
duce the trend seen in the observations. This is a generic 
feature of the cold dark matter paradigm. Any property of 
a galaxy which depends on the mass ratio of the last major 
merger will be scale free. 



4.2 Modifying the Model 

The numerical simulations which motivated the previous 
analysis still leave space for further interpretation. It is e.g. 
not yet fully understood what infiuence gas will have on the 



structure of the merger remnants. Galaxy mergers in nu- 
merical simulations usually start with relaxed spiral galaxies 
even though galaxy formation models predict that mergers 
between ellipt ical galaxies or elliptic al and spiral galaxies 
occur too (Kh ochfar fc Burkertll20o3) . In the following sec- 
tions we therefore study the possible effect of mergers be- 
tween different galaxy morphologies and gas infall on the 
ratio of boxy to disky ellipticals. 



4-2.1 Progenitor Galaxies 

Recent sem i-analytic simu lations of galaxy formation by 
lKhochfar~fc Burkcr^ i2003h showed a remarkable correla- 
tion between the luminosity of elliptical galaxies and the 
morphology of their progenitor galaxies. The most lumi- 
nous elliptical galaxies had preferentially experienced the 
last major mergers between two elliptical galaxies and not 
between two spiral galaxies. Only low luminous galajcies are 
forme d by pure gas rich spiral m ergers. iKormendv fc Bender] 
Jl996tl andEaWeTai] (II993) argued along the same line 
that the most massive elliptical galaxies with boxy isophotes 
are most probably the outcome of a dissipationless merger, 
while less massive ellipticals with disky isophotes result from 
the merger of two gas rich galaxies, which leads to the for- 
mation of a secondary disk in th e remnant. 

The simulations analysed bv lBarned (1199811 : lNaab et alJ 
|l99gD; iNaab fc BurkertJ (|200^ merged spiral galaxies. 
Merge r of spheroids were simulated by iGovcrnato et 
il993r . who showed that the remnant will have most likely 
boxy isophotes. In the following, we will assume that ellipti- 
cals having last major mergers between two bulge dominated 
galaxies {Mtuige O.GMtot) will result in boxy remnants in- 
dependent of the mass ratio. For all the other major mergers 
we adopt the original prescription. 

Figure 13 shows the result of including this dependence 
on the progenitor morphologies for different environments. 
We now find indeed a correlation with luminosity, indepen- 
dent of environment. The most luminous elliptical galaxies 
are mainly boxy and the fraction of boxy ellipticals decreases 
with decreasing luminosity. This trend is an imprint of the 
high fraction of early-type mergers for luminous galaxies. 
iKhochfar fc BurkertI (|20o3) found this fraction to be inde- 
pendent of environment which explains the results in Fig. |3 
However, the gradient which we find is still not as strong as 
observed, although the observational uncertainties are large. 



4-2.2 Disks in Ellipticals 

After the last major merger, elliptical galaxies can continue 
to grow disks; one can identify three main processes of disk 
growth after major mergers. As noted above, gas in tidal 
tails can settle down into a disk after the spheroid has 
formed. We also noted above that we do not model this 
process and instead assume that the gas in the tidal tails is 
accreted to the centre of the remnant at the same time as the 
other present cold gas in the disks of the merging galaxies. 
The second and third process in growing a disk involve the 
cooling of hot gas from halo regions and the accretion of cold 
gas from satellite galaxies in minor mergers, respectively. In 
the following we are going to investigate the influence of the 
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Figure 3. The ratio of boxy to disky ellipticals as function of 
magnitudes as predicted by the simulations. Filled circles show 
the observational data of lBender. Burstein. &: FabeJ 1^993), the 
open symbols show the theoretical predictions for different galaxy 
environments as presented also in Fig. [Tl flower panel). The de- 
pendence on the morphology of the progenitor galaxies in the 
last major merger is included. The mass resolution in all runs is 

Mrmr. = 10^0 Mq . 



Figure 4. Ratio of boxy to disky ellipticals for different mag- 
nitudes found in the simulation. Filled circles show the data of 
iBender. Burstein. &c Faberl ll992l) and open symbols the results 
for the environments presented in Fig. Q (lower panel) including 
the dependence of the isophotal shape on the morphology of the 
progenitor galaxies in the last major merger of ellipticals and the 
formation of a secondary disk by gas infall into the remnant. The 
mass resolution in all runs is Mmin = 10-*^" Mq. 



last two processes mentioned above on the mix of ellipticals 
wit h different isophotal shapes. 

iRix fc Whitel l)l99ff) find that some ellipticals show evi- 
dence for embedded stellar disks with up to 30% of the total 
baryonic mass. It is likely that ellipticals with disks of that 
order will be disky, independent of the details of the last 
major merger. Figure IH shows the results if we include this 
effect in our model. Here we assume that an elliptical will 
be disky if a stellar disk component is present, generated 
by the combined effects of cooling and satellite accretion, 
which contains more than 20% of the total baryonic mass. 
This value is a lower limit based on estimates derived from 
artificially adding an additional stellar disk to a merger rem- 
nant and deriving the isophotal shape of it (Naab, private 
communication) . We ran some tests with different values for 
the lower limit when a elliptical becomes disky, finding that 
for a value of 30% the fraction of boxy to disky elliptical 
did not change. For a value lower than 20% the fraction of 
boxy to disky elliptical galaxies decreases, especially for low- 
luminosity systems which on average have their last major 
merger earlier in contrast to high luminosity elliptical galax- 
ies, which do not have enough time after their last major 
merger to grow a disk that contributes significantly to the 
total light. The fraction of ellipticals that are classified as 
disky because of an additional disk component depends also 
on the environment. In high density environments, the last 
major merger occurs on average at a larger redshift giving 
a new disk more time to grow; this explains the larger frac- 
tion of ellipticals with secondary disk components in high 
density environments on all magnitude scales. Note however 
that these results depend on the efficiency of gas stripping 
in clusters (e.g. Mori & Burkert 2000). 



5 ISOPHOTAL SUB-CLASSES 

The results of the last section suggest that boxy as well as 
disky isophotes have two distinct origins in ellipticals. Each 
formation scenario is expected to leave signatures in the 
structural properties of the remnant. Observatio ns of surface 
brightness profiles also suggest a dichotom y jLauer et alJ 
ll995l:lGebhardt et akllOQ^lFaber et al.ll997^ . Similar to the 
separation into boxy and disky isophotes, ellipticals can be 
classified into core and power-law galaxies, depending on the 
inner slope of the surface brightness profile. High luminous 
elliptical galaxies show shallow inner profiles (core galaxies) 
and low luminous ellipticals show steep inner slopes (power- 
law galaxies). It is still a matter of de bate whether this di- 
chotomy is real or artificial as argued bv lGraham fc GuzmanI 
(2003) who find a continuous linear relation between the 
Sersic index n and the absolute magnitude of the ellipti- 
cal galaxies. For our present analyses we will assume that a 
dich otomy exist . 

iRest et all i200 J) analysed the correlation between 
isophotal shape and the shape of the central surface bright- 
ness profile. They showed that most of the boxy ellipticals 
and disky ellipticals are classified as core and power-law 
galaxies, respectively. Note that not all of the boxy ellipticals 
have a core and not all of the disky ellipticals are power-law 
galaxies, suggesting that each isophotal class might be sub- 
divided into two subclasses, depending on their core prop- 
erties. In the following we investigate the origin of these 
subclasses and test whether our model can account for the 
observation bv lRest et alJ (120011) . 

5.1 Boxy Ellipticals: flat-density cores versus 
power-law cores 

Our model predicts that boxy ellipticals form by two mecha- 
nisms, equal mass mergers of two galaxies and major mergers 
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Figure 5. Ratio of power-law to core boxy elliptical galaxies. 
Stars show the data o f iRest et al1i200lD . The other symbols show 
the theoretical results for the different environments. The mass 
resolution in all runs is Mmin = 10^" Mq. 

{Mgai,i/Mgai,2 < 3.5, Mgai,i ^ Mgaia) of two spheroldally 
dominated galaxies. It is widely accepted that spheroids har- 
bour super - massive black holes fSMBH) which follow a M, — 
a relation llFerrarese fc Merrittl2000l:lGebhardt et al.l200Cl) . 
lMilosavlievl^^^Merri^ ~ il20^ir ' investigated the merger of 
two galaxies harbouring SMBHs and found that the rem- 
nant will have a light profile showing a core in the centre. We 
therefore assume that the major mergers between spheroidal 
dominated ellipticals will lead to the formation of core galax- 
ies with boxy isophotes. Equal mass mergers between disk 
dominated galaxies harbouring no SMBH or SMBH with 
small masses will not lead to the creation of a core due to 
binary black hole merging. In fact, gas infall into the centre 
during the merger and subsequent star formation are likely 
to regenerate cusps. We therefore assume that these mergers 
generate power-law galaxies. 

Fig. 13 shows the ratio of power-law to core boxy ellip- 
tical galaxies at each magnitude for different environments. 
We find no significant environmental dependence. The num- 
ber of core galaxies increases from about 0.1 to 3 times the 
num ber of pow e r-law galaxies from low to high luminosi- 
ties. iRest et al] (1200 j) use three main profile classes: core, 
power-law and intermediate. Following the original work of 
[Pabcr ct al. (1997) we classify in termediate galaxi es as core 
galaxies and compare the data of lRest et al.l (l200lf) with our 
results in Fig. |3 We find that the agreement between the 
simulations and the observations in the overlapping lumi- 
nosity range is very good. Unfortunately, due to the limited 
data, it is not possible to compare the observational data 
with the simulations over a wider range of luminosities. 

5.2 Disky ellipticals: central properties, 

1-component versus 2-component galaxies 

The observations of iRest et al.l (l200ll) show that it is very 
unlikely to find disky ellipticals which are also core galax- 
ies. In fact they find only 1 out of 21. Our simple model 
does not allow for disky ellipticals which have a core. The 
reasoning behind this is that we assume core galaxies form 
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Figure 6. Ratio of the number of disky galaxies formed in major 

mergers with 2 > Mgai,i/Mgai,2 < 3.5 (Mgai,i ^ Mgai,2) to 
disky galaxies which grew a disk with M^isk ^ ^■'^^totMryons 
by gas accretion. The symbols show the results for the different 
environments. No results for magnitudes Mg < —18 in the run 
with Mo = IQl^M© are shown because we did not find any disky 
ellipticals which formed due to gas accretion at these magnitudes. 
The mass resolution in all runs is Mmin = 10-'^'^ Mq. 

by the mergers of two bulge dominated galaxies harbour- 
ing SMBH. This results in the formation of boxy ellipticals. 
In principle, it is possible to generate disky core-ellipticals; 
consider two disk galaxies with significant bulge components 
Mbuig^ O.SMtot. The bulges will harbour SMBH. Unequal 
mass mergers between them would result in disky ellipticals. 
Cores might still form when the massive black holes spiral 
into the centre. At the moment, due to the lack of detailed 
simulations it is not clear whether black hole merging in 
this case would be efficient enough to generate cores. In ad- 
dition, since the combined mass of the progenitor SMBHs 
is smaller than that expected for the SMBH in the remnant 
according to the M, — a relation, the black holes must grow 
significantly in mass by gas accretion in contrast to the case 
of two bulge dominated merging galaxies. Gas infall into 
the central region could trigger central star formation which 
might change the core density distribution, leading again to 
cuspy profiles. Given this complexity, we do not try to model 
the fraction of disky core galaxies, which according to the 
observations is small anyway. 

The fraction of disky elliptical galaxies with massive 
secular disks formed after the last major merger is very sen- 
sitive to both the environment and the amount of residual 
star formation in elliptical galaxies. Recent results from the 
GALEX mission indicate that there is indeed residual star 
formation occuring in elliptical galaxies (S. K. Yi, private 
communication). It is not yet clear where exactly this star 
formation occurs even though there are indications that it is 
in a disk-like substructure. Figure |S| shows the ratio of disky 
ellipticals formed by unequal mass mergers between disk 
galaxies (1-component disky) to disky ellipticals which have 
grown secular massive disks after their last major merger 
by gas accretion (2-component-disky). Most disky ellipti- 
cals are predicted to be 1-component systems that resulted 
from unequal-mass mergers. We find a larger fraction of 2- 
component disky ellipticals in high density environments 
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rather than in low density environments. In addition the 
fraction of 2-component disky ellipticals decreases toward 
higher luminosities. In high density environments, at lumi- 
nosities of Mb ~ —18, the fraction of 2-component disky el- 
lipticals is comparable to the fraction of 1-component disky 
ellipticals. 



6 DISCUSSION AND CONCLUSION 

We used semi-analytical simulations to investigate the for- 
mation of disky and boxy ellipticals. A comparison with ob- 
servations shows that the isophotal shape cannot only de- 
pend on the mass ratio of the last major merger. In this 
case, galaxy properties which depend only on the mass ratio 
of the last major merger should be scale free. The fraction 
of disky-to-boxy ellipticals, however, shows a strong mass 
dependence. To break the mass degeneracy, we propose that 
the isophotal shape should, in addition to the mass ratio, de- 
pend also on the morphology of the progenitors of the last 
major merger. Our model can reproduce well the observed 
trend of the ratio of ellipticals with boxy and disky isophotes 
with luminosity within the observational uncertainties due 
to the low number statistics. In particular, the inclusion of 
the contribution of an additional disk component results in 
a reduction of the fraction of boxy ellipticals mainly at low 
luminosities, while early-type mergers increase the fraction 
of boxy types at high luminosities. The model predicts a 
dependence on the environment once disk formation by gas 
accretion is taken into account. Galaxy assembly by major 
mergers occurs in high density environments at earlier times 
allowing secular disks to grow for longer times and hence be- 
come a more important component in terms of total baryonic 
mass fraction. 

We find two populations of disky ellipticals, those 
formed by unequal mass mergers, called 1-component disky, 
and those who grew new disks and were originally boxy or 1- 
component disky ellipticals, called 2-component disky. The 
fraction of 2-component disky ellipticals is a strong func- 
tion of the luminosity and environment. As massive ellipti- 
cals form late, they do not have time to grow secular disks. 
Disky massive ellipticals therefore are mainly 1-component 
systems, which result from unequal-mass mergers. The disks 
in 2-component disky ellipticals will be on average younger 
than their bulge populations. We predict that such disks will 
be mainly found in cluster environments where, at luminosi- 
ties around Mb ~ —18, half of the disky ellipticals should be 
2-component systems. If these disks would continue to grow, 
they would transform the elliptical into an Sa or Sb galaxy. 
It is not clear whether the galaxy at an intermediate stage 
will look like an SO galaxy or whether SOs are sometimes 
2-component ellipticals seen under certain projections. Ob- 
servations analysing the stellar orbital distribution in disky 
ellipticals will be able to distinguish between the two popu- 
lations of disky ellipticals. 

We predict also two populations of boxy elliptical galax- 
ies, power-law-boxy, made by equal mass mergers of disk 
galaxies, and core- boxy, made by any major merger but be- 
tween two elliptical galaxies. We find that the fraction of 
core-to-power-law galaxies increases toward higher luminosi- 
ties because of a higher fraction of mergers between bulge 
dominated galaxies. However, we do not find an environ- 



mental dependence due to th e self-similar build-up through 
mergers. The observations of lRest et alJ ll200j) agree very 
well with our predictions. A larger observational study will 
be required to test the relation which we find for a larger 
ma gnitude range. 

iBender. Burstein. fc Fabeil (^Q^) argued that ellipti- 
cals can be divided in two main categories, boxy and disky 
ellipticals. Our study shows that to recover observational 
trends we need to make assumptions on the formation 
process of ellipticals which introduce additional sub-classes. 
The boxy and disky sub-classes have different origins which 
are closely related to the hierarchical build-up of structure 
and local gas physics. The two boxy classes differ because 
of the properties of their progenitors which are again 
dependent on the previous merging history. The two disky 
classes differ because of the baryonic gas physics, namely 
the build up of a disk by gas accretion after the last major 
merger. By comparing the predicted fractions of boxy and 
disky sub-classes with observations it will be possible in 
the future to clarify the importance of gas physics coupled 
to star formation and the influence of the global dark 
matter merging history on the origin of early-type galaxy 
morphologies. 
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